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OSCILLATION POWER SPECTRA OF THE SUN AND OF a CEN A: OBSERVATIONS VERSUS MODELS 
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ABSTRACT 



Hydrodynamical, 3D simulations of the outer layers of 
the Sun and a Cen A are used to obtain constraints on the 
properties of turbulent convection in such stars. These 
constraints enable us to compute — on the base of a the- 
oretical model of stochastic excitation — the rate P at 
which p modes are excited by turbulent convection in 
those two stars. Results are then compared with solar 
seismic observations and recent observations of a Cen A. 
For the Sun, a good agreement between observations and 
computed P is obtained. For a Cen A a large discrepancy 
is obtained which origin cannot be yet identified: it can 
either be caused by the present data quality which is not 
sufficient for our purpose or by the way the intrinsic am- 
plitudes and the life-times of the modes are determined 
or finally attributed to our present modelling. Neverthe- 
less, data with higher quality or/and more adapted data 
reductions will likely provide constraints on the p-mode 
excitation mechanism in a Cen A. 

Key words: turbulence, convection, oscillations, excita- 
tion. Sun, a Cen A. 



iButler et"an ('2004') have obtained ultra-high-precision 
velocity measurements of oscillations in a Cen A. From 
those observations Bedding et al. (2004) have derived os- 
cillation amplitudes and also averaged estimates of the 
oscillation life-times. 

The stochastic excitation mechanism has been modeled 
by several authors ('e.g. Goldreich & Keelev 19771 
iGoldreich et al."1994\ 'Samadi & Goupil .2001 ). ' On the 
base of (Samadi & Goupil (200H's theoretical model of 
stochastic excitation and constraints from a 3D simula- 
tion of the Sun, Samadi et al. (2003a) succeeded in mod- 
eling the rates P at which the solar p-modes are excited 
by using a Lorentzian function as a model for the eddy- 
time correlation Xk, instead of using the Gaussian func- 
tion. This result is summarized in Sect.0] 

An open question is then whether or not such non- 
Gaussian model for x, is also appropriate for other solar- 
like oscillating stars and whether the present model of 
stochastic excitation is valid for other stars? This ques- 
tion is addressed on the basis of the recent seismic obser- 
vations of a Cen A. 



1. INTRODUCTION 



2. THE THEORETICAL MODEL OF STOCHASTIC 
EXCITATION 



Solar-like oscillations are stochastically excited by turbu- 
lent convection and damped by several mechanisms. The 
square of the mode amplitude is proportional to the ratio 
between the rate P at which the the mode is excited and 
the mode damping rate -q. The latter is proportional to 
the mode line-width (F) which is inversely proportional 
to the mode life-time (t). 

Providing that measurements of the oscillation ampli- 
tudes and life-time (or line-width) are available it is pos- 
sible to compute P and hence to derive constraints on 
the models of stochastic excitation. Such measurements 
have been avail a ble for the Sun for severa l years (e.g. 
iT.ihhrechri 119881 IChanlin etalJ 119971 HHl). Recently 



The model of stochastic excitation we consider in this 
work is basically that of lSamadi & G oupil' C2001). This 
model provides an expression for the rate P at which a 
given mode with frequency i^o is excited. This expres- 
sion can be written in a schematic form as: 

P(i/o)oc / dm / dvdTi.{SS).i (1) 



The calculation of the rate P{vo) at which a given p-mode 
of frequency is excited, then results from an integra- 
tion over the stellar mass m and local integrations over 
distance r and time T(r) of the mode eigenfunction ^ and 
the correlation product of the excitation sources {SS). 
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The excitation sources, S, are : the turbulent Reynolds 
stress and the advection of the turbulent entropy fluctu- 
ations by the turbulent motions. < SS > is expressed 
in term of the turbulent kinetic energy spectrum E{k, v) 
whe re k is the wav enumber of a given turbulent e lemen t 
fsee ISamadi & Gou t^il 2001). Following Steii^ (dH), 
E{k, v) is split into a spatial component, E{k), and a fre- 
quency component — Xfc(^)' ^Iso designed as the eddy- 
time correlation function — so that 

E(k,y)=E{k)xk{y) ■ (2) 



3 . COMPUTATION OF THE EXCITATION RATES 

P 



We consider a 3D s imulat ion of a Cen A computed with 
IStein & Nordlund ('1998)'s 3D numerical code which 
models only the upper part of the convective zones. The 
grid of mesh points of the simulation is 50 x 50 x 82 large. 
The simulation has a solar metalicity while the real star 
has a metallicity [Fe/H] ~ 0.2. 

We compute the rate P{v) at which p modes are ex- 
cited in a Cen A a ccording to Eq. (1) in the man- 
ner of ISamadi etall ( 12003b.) : The ID stellar model 
consi stent with the simulations of star s is co mputed 
with IChristensen-Dalsgaard & FrandsenI (Il983!) 's stel- 
lar code and assumes the classical mixing-length for- 
mulation of convection. The mixing-length parame- 
ter of the ID model is adjusted in order that the ID 
model matches the simulation. The eigenfunctions (^) 
and their associated frequencies (v) are computed with 
[Christensen-Dalsgaard & Berthomieu (1991)'s adiabatic 
code. The total kinetic energy contained in the turbulent 
kinetic spectrum, E{k), its depth dependence, and its k- 
dependency are obtained directly from the 3D simulation. 
For the eddy time-correlation function Xk, a Gaussian 
function is usually assumed, i.e. it is assumed that two 
dista nt points in the turbulent medium are un-co rrelated 
(e.g. iGoldreich & Keelevlll977t lBalmforfljlT992l) . Here, 
we investigate both a Gaussian and a Lorentzian form for 
Xk- 

For the Sun (see Fig. 1): the rate P{v) at which so- 
lar p modes are excited P{v), was already computed in 
fSamadi et al. (2003b a) on the basis of a 3D simulation of 
the Sun with a higher spatial resolution (253x253x163) 
than that of the present simulation of a Cen A. 



4. EXCITATION OF SOLAR P MODES 



It is shown — on the base of a 3D simulation of the Sun 
— that the Gaussian function usually used for model- 
ing the eddy-time correlation function, Xk, is inadequate 
(Samadi et al. 2003a) . Furthermore the use of the Gaus- 
sian form under-estimates P as shown in Fig. 1 with the 
blue curve. 

On the other hand, a Lorentzian form fits best the fre- 
quency dependence of v as inferred from a 3D simulation 
of the Sun. Computed values of P based on the present 
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Figure 1. Rates P at which solar modes with fre- 
quency n are excited, orange: excitation rates P de- 
duced fivmJh£_seismi£j}bservations the BiSon net- 
work iChaulin et al.il99 S) blue and red: excitation rates 
computed as explained in Sects. 2 and 3. The blue ( red 
resp.) curve corresponds to a calculation in which Xk is 
assumed Gaussian (Lorentzian resp. ). 



model of stochastic excitation and on a Lorentzian func- 
tion for xi- reproduce better the solar seismic observa- 
tions by Chapl in et al. ( 1998) as shown in Fig. 1 with the 
red curve. 

This result then shows that, provided that such non- 
Gaussian model is assumed, the model of stochastic exci- 
tation is — for the Sun — rather satisfactory without ad- 
justment of free parameters in co ntrast with previo us ap- 
proaches dChriste nsen-D alsgaard 1982; B almforth«1992t 
IGoldreich et all 1994t ISamadi et al. 2001) 



5. a CEN A: THE OBSERVATIONS BY BEDDING 
ET AL (2004) 



Solar-like o scillations have re cently been detected in 
a Cen A bv iButler etaP (l2004 with the UCLES spec- 
trograph. The identification of the modes and the de- 
termination o f the mode amplitude s have been recently 
performed bv lBedding et af] (l2004ft . Fig. 2 presents the 
mode amplitudes in Doppler velocity for the 1=0,1,2,3 
modes which have been detected. 

The authors have a lso derived — on the basis of the 
method proposed by Stello et al.'(2004') — an estimation 
of the mode life-times for two different frequency ranges 
(see Table 1). 



Table 1. Life-tim e of the modes as derived by 
\Beddine et al\ for two different frequency ranges. 



Frequency range (/iHz) 1700-2400 2400-3000 
Mode life-time (days) 1.4 1.3 
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a Cen A (Bedding et al 2004) 
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Figure 2. Amplitudes in Doppler velocity of the modes 
detected in a Cen A by Beddins et al. ( 2004). The green 
filled dots correspond to the 1=0 modes, the black ones to 
the 1=1 modes, the red ones to the 1=2 modes and the or- 
ange one to the 1=3 modes. Here the observed amplitudes 
have not been corrected for the mode visibility. 



a CEN A: INFERRING THE INTRINSIC 
AMPLITUDES OF THE MODES FROM 
OBSERVATIONS 



Figure 3. Inferred intrinsic amplitudes of the 1=1 modes 
(black) and of the 1=0 modes (green). The amplitudes of 
the 1=1 modes are corrected for the mode visibility 
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As can been seen in Fig. 2, the measured amplitudes are 
very scattered. Due to visibility effects, intrinsic ampli- 
tudes of 1=3 modes (/3) are expected to be much smaller 
than those of 1=1 (/I) i.e. 13 < II. 

However some observed amplitudes of 1=3 (03) are 
found as large as those of 1=1 (01) but for a given radial 
order, amplitudes for 1=1, 1=3 modes are clearly anticor- 
related as shown in Fig. 2. This suggests that the observed 
amplitude of a given mode 1=3 can be strongly biaised by 
the presence of a neighborhood 1= 1 mode. In contrast, the 
pollution of the 1=1 amplitude by that of a 1=3 one is ex- 
pected to be quite small. Neglecting this pollution effect, 
we sum the amplitudes of 1=3 and 1=1 neighboring modes 
of a given radial order and attribute this amplitude to 1=1 
mode i.e.: Since /3 < II, then 03+01 = /3+/1 ^ II. 
We next correct for the mode visibility. 
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Figure 4. Filled dots correspond to the mode excita- 
tion rates derived from the inferred intrinsic amplitude 
(Fig. 3) of the 1=0 mode (green) and the 1=1 mode (black) 
amplitudes and the estimated mode life-times (Table 1). 
The blue and red curves correspond to the excitation 
rates computed as explained in Sects. 2 and 3. The blue 
(red resp.) curve corresponds to a calculation in which 
Xk is assumed Gaussian (Lorentzian resp.). 



We perform the same with the 1=0 modes: we sum the 
amplitude of 1=0 mode and the 1=2 mode. However this 
determination of the amplitudes of the 1=0 modes must 
be more biased than for the 1=1 modes as the amplitudes 
of the 1=2 modes are not so small compared with those of 
the 1=0 modes. 

Results of this determination of the intrinsic mode ampli- 
tudes are presented in Fig. 3. 



7. a CEN A: COMPARISON OF THE PREDICTED 
EXCITATION RATES P WITH THOSE 
INFERRED FROM THE OBSERVATIONS 



(Table 1) inferred from the observations according to the 
relation: 

P = 2nMrV^ (3) 

where F is the mode line-width related to the mode life- 
time T (Table 1) as F = I/tt/t, M is the mode mass 
calculated from the mode eigenfunctions associated with 
the ID stellar model of a Cen A and finally V is the in- 
trinsic mode amplitude (Fig. 3). 

Theoretical P are calculated according to Eq. (1) and as 
explained in Sects. 2 and 3. 



We derive the 'observed' excitation rates P from the in- 
trinsic mode amplitudes (Fig. 3) and the mode life-times 



Inferred and predicted mode excitation rates of a Cen A 
are presented in Fig. 4. 
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8. CONCLUSION 



For a Cen A as for the Sun, P computed assuming the 
Lorentzian function for the eddy-time correlation func- 
tion Xfe tend to be closer to the observations than those 
assuming the Gaussian function, however the quality of 
the data does not allow to conclude firmly about the best 
function for Xfc. 

In this first tentative comparison, the predictions and the 
observations of a Cen A disagree by a large amount (the 
inferred P are very scattered: the disagreement varies by 
a factor 2 to 7 approximatively). This large discrepancy 
can either be caused by the present data quality which 
is not sufficient for our purpose , by the way the intrin- 
sic amplitudes and the life-times of the modes are deter- 
mined or finally attributed to the present modelling. Data 
of higher quaUty or/and more adapted data reductions are 
really welcome for deriving precise constraints on the p- 
mode excitation in a Cen A. 



